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DIFFUSION  OF  SMALL-SCALE  DENSITY  IRREGULARITIES 
DURING  EQUATORIAL  SPREAD  F 

I.  INTRODUCTION 

Ovar  the  past  few  years  a  considerable  amount  of  effort  has 
been  directed  at  understanding  the  physical  processes  associated 
with  equatorial  spread  F.  Major  advances  have  been  made  both 
experimentally  (Kelley  et  al . .  1976;  Woodman  and  La  Hoz.  1976; 
McClure  et  al .  ,  1977;  Hubo.  et  al .  ,  197$;  Weber  et  al-  ,  1978:- 
Tsunoda  and  Towle,  1979;  Towle ,  1980;  Tsunoda ,  1980, 

Szuszczewicz  et  al . ,  1980;  Tsunoda .  1981;  Rino  et  al. ,  1981; 
Keskinen  et  al.  ,  1981)  and  theoretically  (Haerendel ,  I97U; 

Hudson  and  Kennel,  1975;  Scannapieco  and  Ossakow.  1976;  Chaburvedi 
and  Ossakow,  1977;  Ott ,  1978;  Ossakow  and  Chaturvedi,  1978;  Costa 
and  Kelley,  1978a, b;  Huba  et  al. .  1978;  Ossakow  et  al. ,  1979; 
Keskinen  et  al. ,  1980;  Zalesak  and  Ossakow,  1980;  Sperling  and 
Goldman,  1 9 8 0 ;  Keskinen  et  al .  ,  1981  ;  Huba  and  Ossakow,  1981a ,b ;  for 
a  complete  review  through  1980,  see  Ossakow,  1981)  in  this  area. 
Presently,  much  of  the  research  on  spread  F  is  focussed  on  small- 
scale  irregularities;  that  is  density  fluctuations  occurring  on 
scale  lengths  less  than  several  hundred  meters.  A  brief  overview 
of  this  work  is  as  follows. 

The  first  indication  of  small-scale  density  fluctuations 
being  present  during  equational  spread  F  is  the  3  m  backscatter 
radar  measurements  made  at  Jicamarca  in  the  early  seventies 
( Farley  et  al .  .  1970).  Since  the  mean  ion  gyro-radius  is  'v  U  m, 
these  observations  show  that  turbulence  exists  on  scale  lengths 
smaller  than  the  ion  gyro-radius.  Moreover,  power  spectral  den¬ 
sity  and  radar  backscatter  measurements  (Woodman  and  Basu.  1978) 
Manuscript  submitted  March  26, 1980. 
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suggest  that  the  actual  observed  3  m  radar  backscatter  power  is 
substantially  below  (3-1*  orders  of  magnitude)  the  extrapolated 
in  situ  k“~  spectrum  from  longer  wavelengths  (i.e.,  -^100  m). 

Thus,  the  3  m  irregularities  do  not  appear  to  be  caused  by  a 
turbulent  cascade  of  energy  from  smaller  wave  numbers.  Recently, 
backscatter  measurements  have  been  made  using  the  ALTAIR  ( Hub a 
et  al. ,  1978;  Towle ,  1980;  Tsunoda  and  Towle.  1979;  Tsunoda, 

1981)  and  TRADEX  (Tsunoda .  1980)  radar  systems  at  Kwajalein  and 
have  found  fluctuations  at  1  m,  36  cm,  and  11  cm.  It  appears 
then  that  turbulence  even  exists  at  wavelengths  comparable  to 
the  mean  electron  gyro-radius  (~3  cm). 

In  situ  satellite  (McClure  et  al. .  1977)  and  rocket  (Kelley 
et  al. ,  1976;  Morse  et  al..  1977;  Szuszczewicz  et  al. .  1980)  ex¬ 
periments  have  also  been  performed  during  equatorial  spread  F.  The 
major  result  relevant  to  the  small-scale  irregularities  is  the  ob¬ 
servation  of  steep  plasma  density  gradients  (Costa  and  Kelley, 
1978a;  Kelley  et  al  ,  ,  1981;  Szuszczewicz ,  private  communication, 
1980).  These  sharp  density  gradients  presumably  arise  from  the 
nonlinear  development  of  plasma  macroinstabilities  (e.g., 
Rayleigh-Taylor ) .  Density  gradient  scale  lengths  have  been  ob¬ 
served  as  small  as  30  m,  although  typically  they  are~75  m. 

Also,  correlative  studies  indicate  that  the  strongest  radar  back¬ 
scatter  signals  coincide  with  the  walls  of  plasma  bubbles,  i.e., 
regions  of  steep  density  gradients  (Szuszczewicz  et  al .  .  1980; 
Tsunoda .  1981).  Based  on  these  observations,  the  most  plausible 
explanation  for  the  small-scale  irregularities  is  the  excitation 
of  drift  waves.  In  the  wavelength  regime  £  1  (r^  is  the 

near  ion  gyro-radius),  which  corresponds  to  scale-sizes  greater 
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than  25  m,  the  universal  drift  (collisionless;  Costa  and  Kelley. 
1978a, b),  or  drift  dissipative  ( collisional ;  Goldman  and  Sperling, 
198O;  Huba  and  Ossakov.  1979a)  instabilities  may  be  operative  de¬ 
pending  upon  the  physical  parameters  (e.g.,  collision  freguenc- 
ies).  However,  at  wavelengths  corresponding  to  the  3  m  irregul¬ 
arities,  ion  viscuous  camping  prevents  the  linear  excitation  of 
these  instabilities  for  typical  spread  F  conditions  (Huba  and 
Ossakov ,  1979a).  It  has  been  suggested  that  these  irregulari¬ 

ties  may  be  generated  via  a  parametric  process  driven  by  a  large- 
amplitude,  long  wavelength  mode  (Huba  and  Ossakow,  1979a).  On  the 
other  hand,  the  generation  of  the  1  m,  36  cm  and  11  cm  irregular¬ 
ities  are  probably  due  to  the  lower-hybrid-drift  instability  (or 
possibly  the  drift  cyclotron  instability)  which  exists  in  the  re¬ 
gime  krT  $  1,  where  rT  is  the  mean  electron  gyro-radius  ( Huba 
et  al .  ,  1978;  Huba  and  Ossakow.  1979b;  Huba  and  Ossakow,  1981a, b) 

A  question  naturally  arises:  what  is  the  influence  of  the 
drift  wave  turbulence  on  the  evolution  of  the  plasma?  Laboratory 
experiments  indicate  that  the  dominant  effect  is  anomalous  dif¬ 
fusion  of  plasma  across  the  magnetic  field  to  smooth  out  density 
gradients.  That  is,  particles  interact  with  the  collective  elec¬ 
tric  fields  associated  with  the  instability  and  are  able  to  scat¬ 
ter  across  field  lines.  Thus,  in  general,  drift  instabilities 
act  to  destroy  the  free  energy  source  which  drive  them,  i.e  ,  the 
density  gradients.  Of  course,  the  ionospheric  F  region  is  not 
collisionless  and  the  classical  electron-ion  collision  frequency  can 

be  significant  (i.e.,  v  .  /  0.  ^  1-5  where  v  „  is  the  classical 
0  cfc  i  cl 
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electron-ion  frequency  and  is  the  ion  gyro-frequency)  de¬ 
pending  on  the  value  of  the  density.  Coulomb  collisions  may  al¬ 
so  he  important  in  the  cross-field  diffusion  of  plasma.  Collis¬ 
ions  with  neutral  particles  also  occur  but  are  less  frequent  than 
electron-ion  collisions  at  altitudes  above  300  km. 

Observationally ,  diffusion  appears  to  be  the  dominant 
mechanism  which  smooths  out  the  steep  density  gradients;  those 
such  that  L  ~  several  hundred  meters  where  L  is  the  density 
gradient  scale  length.  By  this  we  mean  the  following.  In  the 
hierarchy  of  small-scale  irregularities  responsible  for  the 
radar  backscatter  measurements,  the  smaller  the  irregularity 


size,  the  steeper  the  density  gradient  scale  length  required  to 
support  it.  Very  sharp  density  gradients  (Ln£  100  m)  are 
necessary  to  excite  the  lower-hybrid-drift  instability  which  is 
responsible  for  the  1  m,  36  cm,  and  11  cm  irregularities.  On 
the  other  hand,  weaker  density  gradients  (Ln>  ICO  m  )  can 
drive  the  longer  wavelengths  modes  (krT.:$  1)  which  are  probably 
necessary  to  generate  the  3  m  irregularities.  Since  the  time 


scale  associated  with  a  diffusion  process  is  t  'v-X^/D 

D 


scale*  length  and  D  is  the  diffusion  coefficient),  the 


( X  is  a 
shortest 


density  gradient  scale  lengths  diffuse  away  first.  Thus,  one 


would  expect  the  smallest  scale  irregularities  to  disappear 


first  in  the  decay  phase  of  equatorial  spread  F  and  this,  in 
fact,  seems  to  be  the  case  ( Baau  et  al .  .  1978;  Basu  et  al .  ,  1980). 

The  purpose  of  this  paper  is  to  examine  both  clsasical  and 
anomalous  diffusion  processes  for  equatorial  spread  F  conditions 


As  such,  we  will  neglect  all  driving  forces,  which  generate 
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equatorial  spread  F  irregularitites ,  and  consider  the  pure  dif¬ 
fusive  decay  of  such  irregularities.  We  focus  on  the  steeper 
density  gradient  scale  lengths  (i.e.,  1  several  hundred 

meters),  since  diffusion  is  much  too  slow  to  be  important  in  the 
evolution  of  the  large-scale  irregularities  (Ln^  1  km).  Our 
main  conclusion  is  that  anomalous  diffusion  is  the  dominant  dif¬ 
fusion  mechanism  for  spread  F. 

The  scheme  of  the  paper  is  as  follows.  In  the  next  section 
we  present  the  diffusion  equation  and  a  discussion  of  the  diffus¬ 
ion  coefficients  considered.  For  the  anomalous  diffusion  coef¬ 
ficient  we  base  our  analysis  on  the  universal  drift  instability. 
In  Section  III  we  discuss  this  choice  of  the  anomalous  diffusion 
coefficient.  Section  IV  contains  the  results  of  a  numerical  a- 
nalysis  of  the  diffusion  equation  for  both  classical  and 
anomalous  diffusion  coefficients.  In  the  final  section  we  apply 
our  results  to  equatorial  spread  F.  We  also  present  similarity 
solutions  to  the  diffusion  equation  in  the  Appendix. 

II  THEORY 

We  consider  the  problem  of  cross-field  diffusion  of  plasma 
in  a  low  g  plasma  (g<<l,  where  g  =  8irn  ( Tg+Ti /B2  )  )  .  The  one¬ 
dimensional  diffusion  equation  can  be  written  as 


where  [Perkins  et  al.,  1973} 


D 


(1) 


(2) 
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\ >  .  is  the  electron-ion  collision  frequency  (either  classical 
ei 

h 

or  anomalous)  and  w  =  (2*irne2/m  P  is  the  electron  plasma  fre- 

P6  6 

quency.  Note  that  we  assume  a  one-dimensional  slab  geometry 

which  is  adequate  for  small-scale  spread  F  irregularities.  The 

one  dimensionality  simplifies  the  analysis,  of  course,  but  it 

is  a  reasonable  assumption  given  that  the  irregularities  during 

ESF  are  generated  by  steep  gradients,  e.g.,  the  walls  of  bubbles. 

Also,  we  assume  T  =  T.  and  have  neglected  electron-neutral 
5  el 

collisions . 

We  choose  the  following  collision  frequencies  for  our 
analysis  : 

c  3/2  _1 

Classical:  v  =v  a(X  /3.5  x  105)(n  /T  )  sec  (3) 

61  0  £  C  6  6 

Anomalous:  v  *v  =  (2ir/9)  |e  r.  .  |  ft  (4) 

61  311  n  JL 1  6 

In  Equations  (3)  and  (1+)  is  the  Coulomb  logarithm 

(X  ■=  23.4-1.15  log  n  +  3.45  log  T  ),  e  =  Id  Bm  n/d  |  •=  1/L 

c  e  °  e  n  1  x1  n 

where  L  is  the  density  gradient  scale  length,  rT.  is  the  mean 

tl  LI 

ion  gyro-radius,  is  the  electron  gyro-frequency,  ng  is  in 
—  3 

cm  and  T  is  in  eV.  The  anomalous  collision  frequency  used 
e 

is  based  upon  the  anomalous  transport  properties  associated  with 
the  universal  drift  instability  ( Gary ,  1980).  We  discuss  this 
choice  shortly.  These  collision  frequencies  lead  to  the  follow¬ 
ing  diffusion  coefficients: 
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Classical : 


D  =  D  .  =  v  .  p2  (n/n  ) 
cl  cl  es  o 


Anomalous  :  D  =  D  =  — 


•n  rLi'P 


■9  n 

es  e 


where  p 


[(Te/Ti)/me]ls/ne  is  an  effective  electron  gyro- 


radius,  n  is  a  normalization  density  and  v  *v  (n  =n  ) . 

o  c  l  cl  e  o' 

From  Eqs.(5)-(7)  it  is  evident  that 


n  v  , 

D  £i 

an  e 


e  rT  . 
n  Li 1 


Thus,  for  typical  spread  F  conditions  it  is  found  that  D  <<D 

cl  an 

and  one  expects  anomalous  diffusion  to  dominate  over  classical 
diffusion . 

We  finally  substitute  Eqs.  ( 5 ) — ( 7 )  into  Eq.  (1)  and  arrive 
at  the  following  diffusion  equations: 


Classical:  -r—  =  v  0  p2  -- 

d  t  Cl  C  S  d 


Anomalous:  — T  =  4"^  8  p2 

3t  9  e  Mes 


(i^  If) 

a^(i£n  rLil  3x) 


We  now  transform  to  dimensionless  variables  and  obtain 


n ,  .  ,  3n  3  ( %  3n  \ 

Classical :  ^  {  ■>  ) 
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(11) 


Clearly,  Eqs  .  (10)  and  (ll)  are  non-linear  partial  differential 

equations.  We  solve  them  numerically  as  initial  value  problems 
in  Section  IV.  They  also  possess  similarity  solutions  which 
offer  some  insight  into  the  scaling  of  diffusion  process.  We 
discuss  these  solutions  in  the  Appendix.  However,  we  now  discuss 
our  choice  of  the  anomalous  diffusion  coefficient. 


III.  ANOMALOUS  DIFFUSION 

As  mentioned  earlier,  the  collective  electric  fields  as¬ 
sociated  with  a  plasma  instability  driven  by  a  density  gradient 
can  cause  particles  to  be  scattered  across  magnetic  field  lines, 
which  leads  to  anomalous  diffusion  of  plasma.  This  process  can 
eventually  smooth  out  the  density  gradient  and  suppress  the 
instability.  Several  estimates  of  the  anomalous  diffusion  co¬ 
efficient  for  this  phenomenon  have  been  given  based  upon  a 
variety  of  physical  arguments.  Perhaps  the  best  known  is 


D  <  y/k  ; 

an  x 


(12) 


which  is  derived  heurist ically  by  Kadomtsev  (1965)  for  a 
strongly  turbulent  plasma  (i.e.,  y  ^  u>r  when  w  =  +  iy). 

More  sophisticated  derivations  of  this  relationshin  have  beer 
presented  although  there  are  criticisms  of  this  estimate. 
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Nevertheless,  Eq.  (l2)  is  used  frequently,  primarily  because  of 
its  simplicity. 

An  alternative  procedure  to  estimate  the  anomalous  dif¬ 
fusion  coefficient  is  based  upon  quasi-linear  theory.  Second- 
order  Vlasov  theory  is  used  to  define  an  anomalous  collision 
frequency  which  can  then  be  determined  from  the  linear  pro¬ 
perties  of  the  mode  and  the  saturation  energy  of  the  instability. 
The  diffusion  coefficient  is  then  found  from  Eq.  (2). 

In  this  paper  we  opt  for  the  second  approach  and  base  our 
analysis  on  the  universal  drift  instability,  i.e.,  collisionless 
drift  instability,  (it  should  be  noted  that  Kelley  et  al . ,  1981, 
have  suggested  that  collisionless  drift  mode  waves  are  the  most 
likely  candidate  for  the  observed  electric  field  and  density 
fluctuation  spectra  in  the  10-100  m  regime).  The  details  of  this 
method  are  outlined  by  Gary  (1980)  and  we  do  not  reproduce  them 
here.  However,  several  comments  should  be  made.  First,  the 
estimate  of  D  given  by  Eq.  (7)  is,  in  fact,  comparable  to 

all 

that  of  Eq.  ( 12) .  Second,  Eq.  (7)  explicitly  contains  the  den¬ 
sity  gradient  which  is  consistent  with  the  notion  that  dif¬ 
fusion  occurs  in  regions  that  are  unstable.  This  dependence  on  L 

n 

leads  to  the  nonlinearity  of  the  diffusion  equation.  Finally, 
other  instabilities,  such  as  the  drift  dissipative  (collisional 
drift  instability)  (Goldman  and  Sperling,  1979)  or  the  lower- 
hybrid-drift  instability,  could  also  be  relevant.  Wnether  the 
universal  drift  instability  is  excited  depends  upon  the  classical 
electron-ion  collision  frequency  v  ,  which,  in  turn,  depends 


9 


upon  the  density.  Typically,  vcj/^i  ^  0*8-5  for  electron  density 
4  6  _  3 

in  the  range  10  -  10  cm  .  For  unstable  modes  in  the  regime 

krT .  <  1  it  is  found  that  the  drift  instability  is  neither 
Li 

purely  collisionless,  nor  purely  collisional  (Huba  and  Ossakov. 
1979a).  Also  estimates  of  the  diffusion  coefficient  associated 
with  the  drift  dissipative  instability  are  comparable  to  that 
of  the  universal  drift  instability  ( Kadomtsev,  1965).  Thus, 
although  collisions  affect  the  universal  drift  instability  for 
typical  spread  F  conditions,  use  of  the  anomalous  diffusion  co¬ 
efficient  associated  with  this  mode  is  Justified  since  it  is 

_1 

qualitatively  accurate  (i.e.,  D“Ln”  )  and  is  quantitatively  ac¬ 
curate  to  within  a  factor  of  2  or  3.  On  the  other  hand,  the 
anomalous  diffusion  coefficient  of  the  lower-hybrid-drift  in¬ 
stability  is  substantially  smaller  than  that  of  the  universal 
drift  instability  and  can  be  neglected  ( Gary,  1980). 

IV.  NUMERICAL  RESULTS 

We  solve  Eqs.  (10)  and  (11)  numerically  as  an  initial  value 
problem  using  a  leap-frog  scheme  for  the  temporal  variation  and 
a  fourth-order,  finite  difference  scheme  for  the  spatial  varia¬ 
tion.  We  choose  as  the  initial  density  profile 

n  (s,  t=»0)  =  (1  -e)  [l  +  e  tanh  s}  (13) 

'Vi 

where  0<e<l,  s  *  x/A,  and  n*n  ( s ) /n  ( s* -®>)  .  Here,  e  determines 
the  magnitude  of  the  density  change  across  the  boundary  layer 
A.  The  boundary  condition  is  3n/3s  *  0  at  the  boundaries 
(i.e.,  j  s  j  >  >  1 )  and  n(|s|  »,t)  =  n(|s|  ®,0).  We  chocse  At 
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sufficiently  small  so  that  the  Courant-Levy  condition  is  sat¬ 
isfied  to  prevent  numerical  instability.  The  grid  size  As  is 
varied  to  insure  that  an  asymptotic  solution  is  obtained  as 
As  •+•  0 .  For  the  results  presented  below  we  choose  e  =  0.9  so 
that  the  density  changes  approximately  an  order  of  magnitude 
across  the  boundary  layer  X.  For  different  values  of  e,  the 
results  remain  qualitatively  the  same  although  quantitative 
changes  occur.  We  also  point  out  that  X  is  a  constant  which 
characterizes  the  initial  width  of  the  boundary  layer.  However, 
X  is  not  the  same  as  the  scale  length  of  the  density  gradient 
Ln .  The  density  gradient  scale  length  is  defined  by 
Ln  =  *[d  C*  n/dsl  and  varies  across  the  boundary  layer.  It 
is  the  density  gradient  scale  length  that  is  critical  to  the 
excitation  of  drift  instabilities  and  not  X. 

A.  Classical  Diffusion 

We  expand  Eq.  (10)  to  obtain 


where  t  *  v  „t  (p  2/X2 ) .  The  solution  to  this  equation  is  shown 
cl  cl  es 

r\j 

in  Figure  1  which  plots  n(s,t)  vs.  s  at  times  t  ^  =  0.0,  0.2, 
and  0.4.  As  is  expected,  the  width  of  the  boundary  layer  is  in¬ 
creasing  with  time.  Note  that  diffusion  is  occurring  more 
rapidly  in  the  region  s>0  than  the  region  s<0.  This  occurs 
because  of  the  density  dependence  of  the  diffusion  coefficient. 
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Since  D  *n  and  the  density  is  greater  for  s>0  than  s<0,  the 
c 

diffusion  coefficient  is  larger  for  s>0  and  faster  diffusion 
takes  place  in  this  region. 

In  figures  2a  and  2b,  we  plot  the  spatial  and  temporal 

variation  of  the  inverse  density  gradient  scale  length.  In 

Figures  2a  we  show  A/L  vs.  s(where  A/L  =  (l/n)  3n/3s)  for 

n  n 

t  *  0.0,  0.2  and  0.4.  The  shortest  density  gradient  scale 

c  x> 

lengths  exist  in  the  region  s<0  due  to  the  (l/n)  dependence  of 
1 

L  ..  However,  it  is  interesting  to  note  that  the  maximum  in- 
n 

verse  gradient  scale  length  actually  increases  during  the  initial 
evolution  of  the  density  profile.  This  is  better  shown  in 
figure  2b,  which  plots  (A/L  )  vs.  t  The  inverse  gradient 

scale  length  increases  ^2056  in  a  time  t  •x*  0.06.  After  this 
time  it  decreases  monotonically . 

The  reason  for  the  initial  steepening  of  the  density  pro¬ 
file  can  be  understood  as  follows.  From  Eq.  (l4)  we  find  that 


n ( s , t+At  )  ^  n(s, t)  +  At 


[ 


(3n/3s) 


+ 


'v 

n 


3  2n/  3  s 


2] 


(15) 


for  A  t  <  <  t .  A  portion  of  the  initial  profile  n(s,T)  is  shown  in 

Figure  3  by  the  solid  line  (not  drawn  to  scale),  and  we  have 

isolated  3  points  (so  ,  si  ,  82).  Note  that  (3n/3s)  >  0  and 

(d2n/ds2)  >  0  for  the  region  shown  so  that  both  terms  in  the 

brackets  of  Eq.  (15)  initially  tend  to  increase  the  density  for 

<\» 

so<s<S2.  At  8=81  the  density  increases  an  amount  Ani  and  at 

f\, 

3=82  the  density  increases  amount  Ani,  in  a  time  At.  However, 
An2>Anj  since  ( 3r./3  s)  >(3n/3s)  ,  and  n  >ni.  This  can  be  shown 

Si  si  * 

easily  if  we  take  n'V/S  in  the  region  So<s<S2. 
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Physically,  since  D«n,  the  diffusion  coefficient  is  larger  at 
s 2  than  si  and  more  particles  diffuse  to  &2  from  s>S2  than  to 


from  s>S] .  Thus,  a  steepening  of  the  density  profile  occurs. 

2 

This  process  continues  until  32n/3s£o  at  s»S2  so  that  Ani  >dn2. 
This  can  be  seen  by  noting  the  evolution  of  n  and  (  A/Ijn)  at 
s  =  s  i  in  Figures  1  and  2a.  At  t  «  C.2  the  density  profile  has 
steepened  (i.e.,  A  /  has  increased  from  its  value  at  x*0.0) 
and  32n/3s2>0.  However,  at  t  -  0.4  the  profile  has  become  less 
steep  (i.e.,  A/L^  has  decreased  from  its  value  at  t  »  0.2)  and 
32n/3s2<  0. 

B.  Anomalous  Diffusion 


We  rewrite  Eq.  (11)  as 


.'v 

3n 

3  T 


an 


(16) 


where  t  -  (W9)(p  2/A2)(r../A)  0  t.  Note  the  explicit  de¬ 
an  es  ui  e 

pendence  on  A/Ln  so  that  diffusion  only  occurs  in  regions  of  the 
plasma  density  gradient  (i.e.,  only  the  regions  which  can  support 
drift  waves).  The  temporal  evolution  of  the  initial  profile 
n(s,T“0)  is  shown  in  Figure  4  for  times  Xfln“0.0,  1*0»  an<*  2.0. 

In  contrast  to  the  classical  diffusion  process,  the  largest 
amount  of  diffusion  is  occurring  in  the  region  s<0.  This  is 
simply  due  to  the  fact  that  Dan  is  largest  in  this  region.  The 
spatial  and  temporal  evolutions  of  (A/L^)  are  shown  in  Figures 
5a  and  5b.  In  Figure  5a,  it  is  clear  that  (A/Ln)  is  more  per¬ 
vasive  for  s<0  than  s>0.  Also  the  maximum  inverse  gradient 


scale  length  is  decreasing  in  time.  This  is  shown  in  Figure 
5b,  which  indicates  that  (A/L  )  decreases  monotonically  in 
time  and  has  decreased  by  <>50f  during  the  time  period  consider¬ 
ed.  This  result  is  consistent  with  the  notion  that  the  effect 
of  the  drift  instability  is  to  smooth  out  the  plasma  density 
gradient . 

V.  APPLICATION  TO  EQUATORIAL  SPREAD  F 

It  has  been  shown  that  both  classical  and  anomalous  dif¬ 
fusion  processes  tend  to  smooth  out  density  gradients.  The 
question  remains  as  to  which  process  dominates  during  equatorial 
spread  F.  The  key  parameter  which  can  answer  this  question  is 
the  time  scale  of  each  process.  In  Table  I  we  contrast  the 
classical  and  anomalous  diffusion  time  scales  for  parameters 
typical  of  equatorial  spread  F.  We  remind  the  reader  that 


Classical :  t 


(a2/p  2)  'A, 

e  s  cl  ct 


Anomalous :  t 


(9/4*)  <A2/pe82>  <*/rL1>  ‘an/ne 


where  (n(s  —  -)),  peg2  -  0^+^ ) /m^2  ,  rLi  -  (T^)^, 

■  eBQ/  n^c  and  A  is  the  initial  scale  length.  We  choose 

5  _  3  + 

T  •  Tl  ■  0.  1  eV,  B  ■  0.3  G,  n(s»-®)-10  cm  and  consider  an  0 
plasma.  For  these  parameters  we  find  that  pgg  *  3.5.  cm, 

o  -1  6  _1 

r.  .  “  4.3  m,  v  *  130  sec  and  Q  ■  5.3  x  10  sec  .  The 
*Li  ’  c£.  e 

time  scales  Eqa.  (17)  and  (18)  are  shown  in  Table  I  for  several 
Initial  scale  lengths.  A*  50  m,  100  m,  200  m,  and  500  m.  The 
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corresponding  minimum  initial  density  gradient  scale  lengths  are 

L  =  m,  80  m,  160  m,  and  U00  m  where  L  *  (  d£/,  jn/dx )”  .  It  should 
n  n 

be  noted  immediately  from  Table  I  that  the  time  scale  for  the  clas¬ 
sical  diffusion  process  is  hours  while  that  for  the  anomalous  dif¬ 
fusion  process  is  minutes.  As  anticipated  from  Eq.  17^,  clas¬ 
sical  diffusion  is  a  much  slower  process  than  anomalous  dif¬ 
fusion  for  the  F  region  ionospheric  plasma.  Moreover,  since 
the  initial  tendency  of  classical  diffusion  is  to  steepen  the 
density  profile,  at  least  an  hour  passes  before  the  density 
gradient  begins  to  weaken  for  the  situation  in  Table  I.  On  the 
other  hand,  the  action  of  the  anomalous  diffusion  process  is  to 
always  smooth  out  the  density  gradient.  For  initially  steep 
gradients  (Ln~100  m)  ,  wave  turbulence  can  double  their  scale 
lengths  in  several  minutes  (t  $  5  min).  However,  several  hours 
are  needed  to  diffuse  density  profiles  to  scale  lengths  on  the 
order  of  a  kilometer.  In  order  to  significantly  diffuse  density 
gradient  scale  lengths  greater  than  a  kilometer,  many  hours  are 

required  (£8  hrs.)  since  t  «X3 

&n  • 

There  is  experimental  evidence  that  suggests  steep  density 
gradients  (Ln^  ®)  relax  in  a  time  <  5  minutes.  Szuszczcwicz 
et  al . ( 1980 )  has  observed  intense  1  m  backscatter  from  the  ALTAIR 
radar  at  Kwajalein  during  the  decay  phase  of  equatorial  spread  F. 

The  most  intense  backscatter  signal  appears  to  decay  away  on  the 
time  scale  ;S  5  minutes.  Presumably,  very  share  density  gradients 
exist  to  produce  the  1  m  density  fluctuations  via  lower -hybrid- drift 
instability  (Ln<  100  m).  Since  anomalous  diffusion  can  smooth 
out  the  Bhar?  density  gradients  on  this  time  scale,  then  it  is 
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■ 


likely  that  the  lower-hybri d-drift  mode  is  also  suppressed  on 

this  time  scale.  This  would  cause  the  1  m  backscatter  signals 

to  weaken  considerably,  consistent  with  observations.  Also,  as 

indicated  earlier,  in  the  decay  phase  of  equatorial  spread  P, 

the  smallest  density  irregularity  scale-sizes  disappear  first, 

( Basu  et  al . ,  1978;  Basu,  1980),  i.e.,  36  cm  backscatter  fades 

before  the  1  m  backscatter,  the  1  m  backscatter  fades  before  the 

3  m  backscatter,  and  so  on,  which  is  also  consistent  with  an 

anomalous  diffusion  process  due  to  drift  waves  since  t  «  A3. 

an 

On  the  other  hand,  large  scale  density  irregularities  (Ln £  1  km) 
are  observed  to  decay  in  several  hours  (Basu  et  al. ,  1980; 

Aarons  et  al. ,  1980)  and  this  process  cannot  be  explained  by 
diffusion,  even  anomalous  diffusion. 

VI.  CONCLUSION 

The  purpose  of  this  paper  is  to  examine  the  diffusion  of 
density  irregularities  during  equatorial  spread  F.  We  have  con¬ 
sidered  both  classical  and  anomalous  diffusion  processes.  A 
major  assumption  in  our  analysis  is  that  any  driving  mechanism 
which  produces  the  density  irregularities  ( S  several  hundred 
meters)  has  ceased.  Presumably  the  small-scale  irregularities 
are  generated  nonlinearly  from  a  macroinstability  (e.g.,  Rayleigh- 
Taylor).  Thus,  we  limit  our  attention  to  the  decay  phase  of 
equatorial  spread  F. 

Experimentally,  diffusion  seems  to  be  the  dominant  process 
which  smoot. .0  out  sharp  density  gradients  (Lr  ^several  hundred 
meters).  This  is  based  on  the  fact  that  the  diffusion  time  is 
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proportional  to  some  power  of  the  diffusion  length  [l’qr. .  (17) 

and  (18)].  Thus,  the  shorter  density  gradient  scale  lengths 
diffuse  away  before  the  longer  ones.  Observationally  this 
is  suggested  by  the  order  of  the  decay  of  radar  backscatter 
signals  (i.e.,  the  shortest  ones  fade  away  first  )  which 
are  due  to  drift  wave  turbulence.  Also  the  time  scale 
associated  with  this  process  is  of  the  order  of  minutes. 

Based  on  our  analysis,  we  find  that  n.11  anomalous  diffusion 
process  is  consistent  with  these  observations  while  classical 
diffusion  is  much  too  slow.  The  anomalous  diffusion  is  based 
upon  drift  wave  turbulence  presumed  to  exist  at  wavelengths 
such  that  krLi~l.  We  have  considered  an  explicit  model  for  the 
anomalous  diffusion  coefficient  based  upon  the  collisionless 
universal  drift  instability.  Although  collisions  probably 
modify  the  dispersive  properties  of  the  mode,  the  anomalous  dif¬ 
fusion  coefficient  is  not  expected  to  be  significantly  different 
(neither  qualitatively  nor  quantitatively).  On  the  other  hand, 
diffusion  of  large  scale  irregularities  (X  ~  1  km)  does  not 
seem  possible  since  the  time  scale  associated  with  this  process 
is  many  hours  (  ,£  8  hrs .  '  .  Another  mechanism  must  occur  to 
smooth  out  these  irregularities,  e.g.,  shorting  out  to  the  E 
region  (which  continuously  builds  up  in  the  morning  hours). 

A  possible  scenario  is  that  anomalous  diffusion  causes  the 
gradient  scale  lengths  Ln  <■>  several  hundred  meters,  wh’ch  drive 
the  radar  backscatter  observed  irregularities  (£  3  m),  to  dis¬ 
appear  on  time  scales  of  the  order  of  minutes.  Radar  backscatter 
observed  irregularities  (  i  3  m),  themselves  without  the  driving 
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gradient  scale  lengths,  can  disappear  by  either  classical  or 
anomalous  diffusion.  Scintillation  causing  irregularities 
(>1  km)  cannot  decay  by  either  classical  or  anomalous  diffusion 
because  it  requires  too  much  time  (>8  hrs.).  However,  these 
long  wavelengths  can  couple  effectively  via  the  high  conductivity 
along  the  geomagnetic  field,  to  lower  regions  (i.e.,  E  regions) 
of  the  ionosphere  (both  north  and  south  of  the  equator).  As 
these  E  regions  build  up  conductivity  they  short  out  the  polar¬ 
ization  electric  field  associated  with  the  long  wavelength  ir¬ 
regularities  and  so  short  out  the  irregularities.  The  time 
scale  for  the  conductivity  buildup  would  be  shorter  than  the 
time  required  for  these  long  wavelengths  to  disappear  by 
diffusion . 
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Fig.  2  —  Spatial  and  temporal  evolution  of  X/L n,  the  inverse  scale  length  of  the 
density  gradient,  for  a  classical  diffusion  process,  (a)  Spat''  1  olution  of  X/L,, 
vs.  8.  (b)  Temporal  evolution  of  maximum  inverse  gradient  length  (X/L„)max. 


TABLE  1.  Classical  and  Anomalous  Diffusion  Time  Scales 
for  Parameters  Typical  of  Equatorial  Spread  F 
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APPENDIX 


We  present  similarity  solutions  to  the  nonlinear  classical 
and  anomalous  diffusion  equations  (Eqs.  (10)  and  (11)). 

1.  Classical  Diffusion 

The  classical  diffusion  equation,  i.e.,  based  on  the 
Coulomb  electron-ion  collision  frequency,  has  received  consider¬ 
able  attention.  A  variety  of  similarity  solutions  exist  depend¬ 
ing  upon  the  initial  conditions  and  the  boundary  conditions. 

One  solution  is  given  in  Longmire  (1963). 

We  consider  the  equation  (which  is  (14)  recast) 


3n 
3  T 


1 

2 


7s7 


and  search  for  solutions  of  the  form 


n 


a  *vi 
t  n 


(s/Tb) 


(Al) 


(A2) 


where  a  and  b  are  constants,  Defining  a  new  variable 


£  =  s/  T 


(A3) 


and  substituting  Eqs.  (A2)  and  (A3)  into  Eq.  (Al)  yields 


% 

a- 1  ,  a- 1  3n 

at  n  (  C  )  -bt 


i  2a'2b  a** 


2  T 


(A4) 


We  require 


a  ■  2b  -1 


(A5) 
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so  Chat  the  powers  of  t  cancel.  In  order  to  determine  a  and  b 
we  Impose  the  added  condition  that  the  plasma  can  diffuse  freely 
(i.e.,  Tods  *  constant).  This  leads  to  the  requirement  that 


a  +  b  =  0 

Thus,  from  Eqs.  (A5)  and  (A6)  we  obtain 

a  --  1/3;  b  =  1/3 
and  Eq.  ( A4 )  becomes 


l 

3 


% 

n  (O 


1  3n  ,  I  3  2  n2 
3  4  3£  ”  2  IT7 


(A6) 


(A7) 


.(A8) 


The  solutions  of  Eq .  (A8)  require  n  vanish  at  some  point  so 

that  the  plasma  occupies  a  finite  region  at  any  time.  We  choose 
s  -  ±  sq  (the  vanishing  point)  at  t  *  tq.  The  density  is  then 


given  by 


(A9) 


the  initial  scale  length  is  s 


the  scale  length  to  double  is 

AT 


cl 


'  I 


1  and 


1  (x'vA),  the  time  it  takes  for 


(A10) 


or 


•*,.  '  I- 


es 


c  £ 


(All) 
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which  is  comparable  to  the  numerical  results.  It  is  important 
to  recognize  a  major  difference  between  the  similarity  solution 
Eq.  (A9)  and  the  ionospheric  problem.  Equation  (A9)  considers 
the  diffusion  of  a  density  enhancement  as  opposed  to  a  density 
depletion  which  occurs  in  the  ionosphere.  One  difference  is 
the  initial  steepening  of  the  density  profile  for  a  depletion. 
This  effect  does  not  occur  for  an  enhancement.  However,  the 
asymptotic  solutions  of  both  problems  yield  comparable  time 
scales . 

2.  Anomalous  Diffusion 

We  consider  the  equation 


( A 1  2  ) 


and  look  for  solutions  of  the  form 


n  -  1  -  5  exp  (-<Kx ) )  (A13) 

o  b 

where  x  ■  s  /t  and  0<5<1  substituting  Eq  (A13)  into  Eq .  (A12), 

we  find  a  -  3  and  b  *  I  where  *Hx)  satisfies  the  nonlinear  dif¬ 
ferential  equation 


X 


4  d^ 
3  dX 


2  X 


d2H> 

dX2 


( A 1  4  ) 


For  large  x»  the  solution  to  Eq.  (A14)  is 


=  X 


»s 


(A15) 
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or 


n 


1  -<$  exp  [  -  (sS/t)*5  ] 


(A16) 


This  solution  more  closely  resembles  the  ionospheric  problem 

since  it  represents  a  density  depletion  being  filled  with  plasma. 

Equation  (A16)  predicts  that  the  scale  length  of  the  density 

1/3 

sheath  broadens  at  a  rate  proportional  to  t  .  This  is 

b  n 

slightly  faster  then  the  results  shown  in  Figure  5  for  (A/L  ) 

n  max 

vs.  x  which  is  not  surprising  owing  to  the  approximations 
an 

made  in  obtaining  Eq.  (A16). 
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ATTN  RD/SI,  RM  5648,  HQ  BLDG 
WASHINGTON,  D.C.  20505 

Ola  Attn  OSI/PSID  Rm  5F  19 

DE°AfiTMENT  OF  COMMERCE 

National  Bureau  of  St  awards 

WASHINGTON,  D.C.  20234 

(All  Corres:  attn  Sec  Officer  for) 

OIcy  attn  R.  Moore 

Institute  for  Telecom  Sciences 
National  Telecomminications  8  iNFo  itoMiN 
BOULDER,  Co  80303 

Ola  ATTN  A.  JEAN  UNCLASS  ONLY) 

Ola  ATTN  W.  JTLAUT 
Ola  ATTN  D.  Crombie 
Ula  ATTN  L.  BERRY 
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NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
Environmental  Research  uboratories 

DEPARTMENT  OF  COMWRCE 

BOULDER.  CO  803u2 

OIcy  mttn  R.  Grubb 

OICY  ATTN  AERONOMY  LAB  G.  KEIO 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

aerospace  Corporation 
P,  0.  Box  92957 
los  Angeles,  Ca  9C009 

Ola  Attn  I.  Garflnxel 
Ola  Attn  T.  Salmi 
Ola  Attn  v.  josephson 
Ola  Attn  s.  Boner 
Ola  ATTN  N.  STOCXWELL 
OIcy  Attn  d.  (Xsen 
Ola  ATTN  Slf  A  FOR  PWW 

analytical  Systems  Engineering  Corp 
5  Old  Concord  Road 
BURLINGTON,  MA  01803 

Ola  attn  radio  Sciences 

Berkeley  research  associates,  Inc. 
p.  0.  box  983 
BERKELEY,  CA  94701 

Ola  ATTN  J.  WORKMAN 

BOEING  COWANY,  The 
P.  0.  Box  3707 
Seattle,  wa  98124 

Ola  attn  o.  Keister 

Ola  ATTN  U.  MURRAY 

Ola  attn  g.  hall 
Ola  attn  j.  Kenney 

California  at  San  Diego,  Univ  of 

P.O.  aox  6049 

San  Diego,  CA  92106 

Brown  Engiwering  Cowany,  Inc. 
Cummings  Research  Park 
Huntsville,  al  35807 

Ola  attn  Romeo  a.  Deliberis 

Charles  Stark  Draper  Laboratory,  Inc. 
555  Technology  Square 
Cambridge-  Ma  02139 
Ola  mttn  0.  o.  Cox 
OIcy  attn  j.  P,  Gilmore 


COMSAT  laboratories 
linthicum  Road 
Clarksbirg,  Md  20734 
Ula  ATTN  u.  HYDE 

Cornell  university 

Department  of  Electrical  Engineering 
Ithaca.  NT  14850 

Ola  mttn  0.  T.  Farley  jr 

CLECTROSPACE  SYSTEMS,  INC. 
box  1359 

RlGtARDSON,  TX  75080 

OIcy  attn  n.  logston 

Ola  attn  Security  (Paul  Phillips) 


General  Electric  Cowan/ 

Space  Division 
valley  Forge  Space  Center 

GODDARD  BLVD  KING  OF  PRUSSIA 
P.  0.  Box  8555 
Philadelphia,  Pa  19101 

Ola  Attn  M.  h.  Bortner  Space  Sci  Lab 

General  Electric  Cowany 
P.  0.  Box  1122 
Syracuse-  «Y  13201 

Ola  mTTnF.  RE! BERT 

General  Electric  Cowaw 
Tewo-Center  for  Advanced  Studies 
516  State  Street  (P.O.  Drawer  ou) 

Santa  Barbara.  Ca  93102 
Ola  attn  DASIAC 
Ola  Attn  Don  Chajcler 
01a  mttn  Tom  barrett 
ula  mttn  Tim  Stephans 
ula  mttn  warren  S.  nnapp 

Ola  ATTN  WILLIAM  IXHAMARA 

Ola  attn  B.  Gambill 
Ola  MTTN  MACK  STANTON 

ueneral  Electric  Tech  Services  Co.,  Inc. 
hMES 

Court  Street 
Syracuse,  NY  13201 

Ola  mttn  G.  millman 

GENERAL  RESEARCH  CORPORATION 
Santa  Barbara  Division 
P.  0.  box  6770 
Santa  Barbara,  Ca  93111 
Ola  ATTN  JOHN  ISE  JR 
Ola  ATTN  JOEL  GARBARINO 

Geophysical  Institute 

UNIVERSITY  OF  ALASKA 

Fairbaaks,  ak  99701 

(All  Class  attn.-  security  officer) 
Ola  Attn  T.  N.  Davis  (Uncl  Only) 

Ola  attn  Neal  Brown  (Uncl  Only) 

Ola  attn  Technical  Library 

GTE  Sylvania,  Inc. 

Electronics  Systems  GRP -Eastern  Div 
77  A  Street 
Needham,  ;ia  02194 

OIcy  attn  Marshal  Cross 

Illinois,  university  of 
107  Coble  hall 
150  Davewort  House 
Chawaign,  IL  61820 

(All  Corres  attn  Dan  McClella®) 

Ola  For  K.  yeh 

Institute  for  Defense  Analyses 
400  Army-navy  Drive 
ARLINGTON,  VA  22202 

Ola  attn  j.  m.  aein 
Ola  attn  Ernest  Bauer 
Ola  attn  Hans  wolfhard 

Ola  ATTN  JOEL  bENGSTON 


ESL  Inc. 

m95  JAVA  JRIVE 

Sunnyvale,  Ca  9408b 

Ola  attn  j.  Roberts 
OIcy  mttn  james  Marshall 
tj  iCY  MTTN  C.  W.  ArETTIE 
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hSS.  Inc. 

2  Alfred  Circle 
BEDFORD.  ma  01730 

Ola  Attn  Donald  hansen 

Intl  Tel  &  Telegraph  Corporation 
500  Washington  Avenue 
Nutley,  fU  07110 

Ola  Attn  Technical  Library 


Ratc  Corporation,  The 
1700  Main  Street 
Santa  Monica,  Ca  50*06 
Ola  Attn  Cullen  Crain 
Ola  Attn  Ed  Bedrozian 

Riverside  Research  Institute 
80  nest  End  avenue, 
new  York,  nY  10023 


JAYCOR 

1*01  Cahino  Del  Mar 
Del  Mar.  Ca  92014 

Ola  Attn  S.  R.  uoldhan 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
jcans  hopkins  Road 
Laurel.  Md  20810 

Ola  attn  Document  Librarian 
Ola  Attn  Thomas  Potemta 
Ola  Attn  jow  Dassoulas 

locked  Missiles  8  Space  Co  Inc 
P.  0.  Box  504 
Sunnyvale,  Ca  94088 

Ola  Attn  Dept  oQ-12 
Ola  Attn  D.  R.  Churchill 

LOCKHEED  MISSILES  AIC  SPACE  CO  INC 
3251  hancwer  Street 
Palo  Alto,  U  94304 

Ola  Attn  Martin  Walt  Dept  32-10  _ 

Ola  ATTN  RICHARD  G.  JOHNSON  DEPT  32"12 
Ola  Attn  w.  l.  Imhof  Dept  32-12 


Science  applications.  Inc. 

P.  U.  BOX  2351  ^ 

La  Jolla,  Ca  92038 

Ola  ATTN  LEWIS  M.  LINSON 
Ola  ATTN  UANIEL  a.  HAH. IN 
Ola  ATTN  D.  SACHS 
Ola  ATTN  t.  A.  STRAKER 
Ola  Attn  Curtis  a.  Smith 

nlrv  ATTN  iacx  McDougall 


Raytheon  Co. 

528  Boston  Post  Road 
Sudbury,  Ma  01776 

Ola  Attn  Barbara  Adams 

Science  applications.  Inc. 

1710  Goocridge  Dr. 

MCLEAN,  VA  22102 

Attn:  J.  Cockayne 

Locxheed  Missile  &  Space  Co.^Inc. 
Huntsville  Research  8  £ngr.  Ctr. 
4800  Bradford  Drive 
Huntsville,  Alabama  35807 
attn:  dale  H.  Davis 


KAMAN  SCIENCES  CoRP 
P.  0.  BOX  74b3 
Colorado  springs,  lo  a0933 
Ola  ATTN  1 .  MEAGHER 


Liwcabit  Corp 
10453  ROSELLE 
San  Diego.  Ca  92121 

Ola  Attn  Irwin  jacobs 


McDonnell  Douglas  Corporation 
5301  bolsa  avenue  _ 

HUNTINGTON  BEACH.  CA  92647 

Ola  attn  n.  Harris 

OIcy  attn  J.  Houle 

Ola  attn  ueorge  mroz 

ula  attn  w.  ulson 

Ola  attn  R.  w.  halprin 

r,Vv  JTTW  TcrviMirii  I  1RRARY  Services 


M.I.T.  LINCOLN  LABORATORY 

P.  0.  BOX  73 
LEXINGTON.  MA  02V3 

Ola  attn  David  m.  Towle 

Ola  nTTN  P.  WALDRON 
Ola  Attn  l.  LOUGH,  in 
Ola  attn  D.  Clark 

hartin  Marietta  Corp 
URLANOO  DIVISION 
P.  j.  Box  5837 
JRLANDO.  ?L  32805 

Ola  ATTN  R.  HEFFNER 

Physical  dynamics  Inc. 

P.  0.  Box  3027 
Bellevue,  wA  98009 

Ola  attn  J.  fremouw 

Physical  uynamics  Inc. 

P.  j.  nox  10367 
OAKLAND.  Ca.  94610 
attn:  a.  Thomson 


Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  ua  j3101 
Ola  attn  P.  riscHtR 
Ola  Attn  w.  P.  Crevier 
OIcy  attn  Steven  L.  uutsche 
01a  attn  D.  Sapped i eld 
Ola  Attn  R.  dogusch 
Ola  Attn  R.  hendrick 
Ola  attn  Ralph  Kilb 
Ola  Attn  Dave  SotcE 
OIcy  Attn  P .  Pajen 
OIcy  Attn  M.  Scheibe 
01a  attn  Conrad  l.  longmire 
Ola  Attn  warren  a.  Schjjeter 


Mitre  Corporation.  The 
P.  0.  Box  208 
Bedford,  ma  o1730 

Ola  Attn  john  Morganstern 

Ola  ATTN  G.  HARDING 
OIcy  Attn  C.  t,  Callahan 


r  i  d  associates 
P.  0.  sox  9695 
marina  Dei  Rey,  Ca  90291 

OIcy  attn  Forrest  Gilmore 
ola  Attn  bryan  Gabbard 
OIcy  Attn  william  b.  vuight  Jr 
Ola  attn  Robert  f.  lecevier 
Ola  attn  william  j.  Karzas 
ola  attn  h.  Ory 
OIcy  attn  C.  MacDonald 
Ola  attn  R.  Turco 


mitre  Corp 

WESTGATE  RESEARCH  PARK 
1820  Dolly  .Madison  Blvd 
MCLEAN,  VA  22101 

Ola  ATTN  W.  HALL 

Ola  attn  w.  poster 


ACIFIC-SIERRA  RESEARCH  CORP 
456  Cloverfiecd  Blvd. 
ant a  Monica,  Ca  90404 

Ola  attn  6.  C.  Pi eco  JR 


39 


Pennsylvania  State  University 
Ionosphere  Research  lab  t 
318  Electrical  Ensincerinb  east 
University  Park.  PA  16802 

(NO  CLASSIFIED  TO  THIS  ACCRESS) 

oicy  Attn  Ionospheric  Research  lab 


photonetrics.  Inc. 

4ag  Harriett  Road 
LEXINGTON,  ru  02173 

OIcy  Attn  Irvins  L.  kofsky 


Technology  International  Coap 
75  Wiggins  Avenue 
Bedford.  ha  01730 

OIcy  Attn  w.  P.  Boouist 


TRW  Defense  &  Space  Sys  Group 
One  Space  Park 
kedondo  Beach.  Ca  S0278 
01 CY  ATTN  «.  K.  PLEBUCH 
OIcy  Attn  S.  Altschuler 
01cy  attn  D.  Dee 


SRI  International 
333  ravenswooo  Avejj* 

NEN.0  Park,  la  9RQ2S 

OIcy  Attn  Donald  neilson 
Ola  Attn  Alan  Burns 
Ola  attn  G.  Smith 
Ola  attn  l.  L.  Cobb 
Ola  Attn  David  A.  Johnson 
Ola  Attn  Walt®  G.  Chesnut 
Ola  attn  Charles  l.  Rino 
Ola  Attn  Walt®  Jaye 
Ola  attn  m.  Baron 
Ola  attn  Ray  L.  uadabraao 
Ola  attn  G.  Carpenter 
Ola  attn  g.  price 
Ola  attn  j.  Peterson 
Ola  attn  R.  h*e.  JR. 

Ola  attn  v.  Gonzales 
Ola  attn  D.  McDaniel 
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IONOSPHFRIC  MTORING  DISTRIBUTION  LIST 
IKCLASSIPIED  ONLY 

PLEASE  niSTRIAUTE  ONE  COPY  TD  EACH  OF  THE  FOLLOWING  PEOPLE: 


Advanced  Rf search  Projects  Agency 

(ARPA) 

Strategic  Technology  Office 
Arlington,  Vipginia 

Capt.  Donald  m.  Levine 


Naval  Reseapch  Laboratory 
Washington,  D.C.  20375 

Dr,  P.  Mange 
Dr.  R.  Meier 

Dp.  5.  Szuszczewicz  -  Code  4127 
or.  J.  Goodman  -  Code  7560 


Science  applications,  Inc. 
1250  Prospect  Plaza 
La  Jolla,  California  92037 

Dr.  D.  A.  Hamlin 
Dp.  L.  Linson 
Dr.  D.  Sachs 


Director  of  Space  and  Environmental 
Laboratory  NCAA 
Boulder,  Colorado  8030? 

Dr.  A.  Glenn  Jean 
Dr.  G.  w.  Adams 
Dr.  D.  n.  Anderson 
Dr.  K.  Davies 
Dr.  R.  F.  Donnelly 


A.  F.  Geophysics  Laboratory 
L.  G.  hansom  Field 
Bedford,  Mass.  01730 

Dr.  T.  Elkins 
Dr.  w.  Swider 
Mrs.  R.  Sagalyn 
dr.  J.  m.  Foppes 
Dp.  T.  J.  Keneshea 
hr.  j.  Aarons 


Office  of  Naval  Research 
ROO  North  Quincy  Street 
Arlington,  Virginia  22217 

Dr.  H.  Mijllaney 


Commander 

naval  Ocean  Systems  Center 
San  Diego.  California  92152 

M>.  R.  Rose  -  Code  5321 


U.  S.  Arm/  Aberdeen  Research  and 
Development  Center 
Ballistic  Research  Laboratory 
Aberdeen.  Maryland 

Dr.  J.  Heimerl 


Commander 

Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  D.C.  20360 

IP.  T.  Czuba 


Harvard  University 
Harvapd  Square 
Cambridge,  Mass.  02138 

Dr.  M.  B.  Mcelroy 
Dr.  R.  Lindzen 


Pennsylvania  State  University 
University  Park,  Pennsylvania  16802 

Dr.  J.  S.  Nisbet 
Dr.  p.  R.  Rohrbaugh 
Dr.  L.  a.  Carpenter 
Dr.  m.  lee 
Dr.  p.  Djvany 
Or.  p.  Bennett 
Dr.  f,  klevans 


